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We perform numerical studies of the wave packet propagation through open quantum billiards 
whose classical counterparts exhibit regular and chaotic dynamics. We show that for t < th (th 
being the Heisenberg time), the features in the transmitted and reflected currents are directly 
related to specific classical trajectories connecting the billiard leads. In contrast, the long-time 
asymptotics of the wave packet dynamics is qualitatively different for classical and quantum billiards. 
In particularly, the decay of the quantum system obeys a power law that depends on the number of 
decay channels, and is not sensitive to the nature of classical dynamics (chaotic or regular). 
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Low-dimensional nanometer-scaled semiconductor 
structures, quantum dots (sometimes called the quan- 
tum electron billiards) represent artificial man-made 
systems which are well suited to study different aspects 
of quantum- mechanical scattering jjj. A majority of 
studies of electron transport in such systems have been 
mainly focused on the stationary electron dynamics. 
In recent years, however, interest in temporal aspect 
of quantum scattering has been renewed This 
includes e.g., studies of the time delay distributions || 
and correlation decay in quantum billiards and related 
systems Furthermore, many core starting points in 
the description of the stationary scattering rely heavily 
on the properties of the system in the time domain. 
In particular, the semiclassical approach exploits the 
difference between the classical escape rate from the 
cavities with chaotic and regular (or mixed) dynamics 
(exponentially fast e -7 * for the former vs. power-law 
for the later Q). This difference in the classical 
dynamics translates into the the difference in observed 
transport properties (statistics of the fluctuations 
H 0, ||, a shape of the weak localization ||, etc.). 

On the other hand, the quantum mechanical (QM) ap- 
proaches predict qualitatively different, universal power- 
law escape rate from the cavity [JO] , 

dP(t)/dt ~ r^" 1 , (1) 

where P(t) is the survival probability, M is the number of 
decay channels and j3 = 1 (2) for the system with (with- 
out) time-reversal invariance p2| . The non-exponential 
decay of a quantum system with chaotic classical dynam- 
ics has been indirectly demonstrated Jl2| in a microwave 
stadium billiard. 

The QM power-law delay time for the chaotic cavity is 
expected to deviate from the semi-classical (SC) decay at 
times of the order of the Heisenberg time th = h/A, A 
being the mean level spacing of the cavity flc[ ]. At the 
same time, the difference in the classical decay of the 



chaotic and regular/mixed cavities often becomes dis- 
cernible only after many bounces at the times which 
often exceed th Nevertheless, the SC predictions 

are widely used in experiment to distinguish between the 
chaotic and regular/mixed dynamics in quantum billiards 
0. Is it thus possible to reconcile the SC and QM ap- 
proaches, or should some of the SC predictions be used 
with certain caution or even be revised? Does the long- 
time decay asymptotics of the quantum systems depend 
on the underlying classical dynamics (chaotic or regu- 
lar)? Motivated by these questions we, in this paper, per- 
form direct quantum mechanical calculations of the pas- 
sage of electron wave packets through two-dimensional 
electron billiards. 

To the best of our knowledge, all of the studies of wave 
packets dynamics in open systems presented so far, have 
been mostly restricted to (a) quantum limit where a char- 
acteristic size of the system L was of the order of the 
average wavelength of the wave packet (A) and (b) to an 
initial stage of the wave packet evolution t < 1 (where t 
is in units of the traversal time) . The time-dependent so- 
lution of the Schrodinger equation was typically obtained 
on the basis of direct schemes approximating the expo- 
nential time propagator |p^| . With such methods the 
task of tracing the long-time evolution of a wave packet 
in a realistic quantum dot would be forbiddingly expen- 
sive in terms of both computing power and memory. In 
the present paper we thus implement a spectral method 
based on the Green function technique [HJ , which allows 
us (a) to reach a semi-classical regime (A) 3> L and (b) to 
approach a long-time asymptotics t ^> 1 corresponding 
to 10 4 — 10 5 bounces of a classical particle in a billiard. 
We found that during the initial phase t < th (which 
in our case corresponds to ~ 10 — 20 classical bounces 
in a billiard), the QM decay closely follows the classi- 
cal one, such that all the features in the QM current 
leaking out of the billiard can be explained in terms of 
geometry-specific classical trajectories between the leads. 
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FIG. 1: A square and a quarter-stadium shaped billiard con- 
nected to semi-infinite leads; L = R, L/w = 8. The half-width 
of the wave packet a — OAL; at t = the wave packet is dis- 
tinct from zero in the interval of 3L. The average wavelength 
of the wave packet (A) = 2n/{k) = 0.8w. 



When t > tjj, the QM dynamics starts to deviate from 
the classical one, with the decay rate obeying a power 
law that depends on the number of decay channels only, 
irrespective of the nature of classical dynamics (chaotic 
or regular). We thus conclude that quantum mechanics 
smears out the difference between classically chaotic and 
regular motion. 

We have studied the temporal evolution of wave pack- 
ets in square, Sinai, and stadium billiards of various 
shapes. All of them exhibit similar features and we 
thus present here the results for two representative ge- 
ometries, a square (which is classically regular) and a 
quarter-stadium (which is classically chaotic), see Fig. 
[j]. The billiards arc connected to two semi-infinite leads 
that can support one or more propagating modes. Mag- 
netic field is restricted to zero. We assume a hard wall 
confinement both in the leads and in the interior of bil- 
liards. Dynamics of the wave packet is governed by the 
time-dependent Schrodinger equation 



fo^-H) liK*)> = o, 



(2) 



where H is the Hamiltonian operator and \ip{t)) is the 
wave function. To study the time evolution of the initial 
state we follow St0vneng and Hauge (l5) and perform 
the Laplace transform of Eq. (g) followed by the integra- 
tion by parts. Changing variables in the Mellin inversion 
integral we obtain 



\m) = 7T 



oo+iO 



dzG(z)\ip(0))e 



— izt/h 



(3) 
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where we have introduced the Green function operator 
G(z) = (z — H)^ 1 and taken into account that all the 
poles of the Green function lie in the lower z-plane. With 
the help of Eq. (^), the calculation of the temporal evolu- 
tion of the initial state is effectively reduced to the com- 
putation of the Green function of the Hamiltonian op- 
erator H in the energy domain. To compute the Green 
function we discretize the system under consideration, in- 
troduce a standard tight-binding Hamiltonian and make 
use of the modified recursive Green-function technique 
described in details in |fL6|. 



Let us consider a minimum-uncertainty wave packet of 
the average energy E which enters a billiard from the left 
lead in one of the transverse modes a. We thus write the 
initial state in the left lead at t = in the form 



|Va(0)) = £CXI 



(2tt)VV? 



(4) 



(5) 



y E/u is the average wave vector 
where u — h 2 /2m*a 2 and m 

\2 7,a2 , 7„q 2 



where w is the width of the leads (measured in units 
of a lattice constant a), /" = y/2/wsm^an/w) is the 
eigenfunction of the transverse motion; (k) = 2ir/(\) = 

(in units of a -1 ), 
being the effective mass; 
(k) 2 = kf z + fcj 2 , kf = na/w, with fcj? and kf being 
the longitudinal and transverse wave vectors respectively. 
The matrix element (m, n\i))) defines the probability am- 
plitude to find the electron on the site (m,n). After the 
wave packet enters the billiard, it will leak out through 
both of the leads in all the available modes j3. The wave 
function in e.g. the right lead can then be written in the 
form 



|tfa(*)> N?>. 
in ■> 



(6) 



n •> a 



f a gives a probability to 



where \rr, 

find a particle on the slice m in the transverse mode /?, 
provided the initial state enters the billiard in the mode 
a. Discretizing a standard expression for the quantum- 
mechanical current, j(x,y) — ifi/2m* (ifiVip* — ip*Vip) 
we, using Eq. ([|), obtain the following expression for the 
total current J = J dyj(x,y) through the slice m in the 
leads expressed via coefficients c™ Q (t) 
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(7) 



Note that the quantum-mechanical current is related to 
the survival probability in the billiard Eq. (Q) by the 
obvious relation 



dP(t)/dt = -Ji{t) - J r (t), 



(8) 



where Ji(t) and J r (t) stand for the currents flowing into 
the left and right leads respectively. Note that the func- 
tion Ji (t) + J r (t) has a meaning of the distribution of the 
time delays in the billiard [|| We calculate coefficients 
c™^ 1 in Eq. (Q) by computing a matrix element (mf3\ip) 
using Eqs. @-@ 
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dzJ2G^'(z 



-izt/h 



(9) 



where G™™ (z) stands for the matrix element 
(mj3\G(z)\m'a) of the Green function of the whole 
system (billiard and semi- infinite leads). 
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FIG. 2: Quantum mechanical reflected and transmitted cur- 
rent through a square billiard (upper and lower solid curves 
respectively). Dashed lines indicate corresponding classical 
currents. The wave packet enters the billiard in the second 
mode a = 2; (k) — 2.5w/n. The insets show classical reflected 
and transmitted trajectories; the numbers in the parenthe- 
sis are the so-called winding numbers indicating how many 
times an electron traverses the billiard in the longitudinal 
and transverse directions. The inset to the right show an an- 
gular distribution of injected electrons P{6) calculated in the 
Fraunhofer approximation for a — 2. The Heisenberg time 
th is indicated by a dot-dashed line. The curves are shifted 
for clarity. 



A special care has been taken to ensure the reliabil- 
ity of the results of the numerical simulations. This, in 
particularly, includes the thorough control of the conser- 
vation of the total current. We have also calculated the 
temporal evolution of the wave packet in the infinite lead 
of the width w and found an excellent agreement with 
the analytical results for a ID lattice jlq| . Finally, the 
developed method reproduces correctly the conductance 
quantization of the quantum point contacts. 

Let us first concentrate on the initial phase of the 
wave packet dynamics. Figure || shows the quantum- 
mechanical current for the square billiard in the time in- 
terval t < 15 (here and hereafter we measure time in units 
of the traversal time t tr = L/(v), (v) = L/(h(k)/m*); 
the current is measured in units of h/2m*a). Our choice 
of the parameter of the wave packet er and (k) ensures 
that the spreading of the wave packet becomes noticeable 
only after relatively long time t > 50. The initial period 
of time < t < 3, when the current through the left lead 
Ji(t) is negative, corresponds to a buildup phase when 
the wave packet enters the billiard. Having entered the 
billiard, the wave packet starts to leak out through the 
leads and the calculated QM currents Ji , J r show a series 
of pronounced peaks. To outline the origin of these peaks 
we calculate the leakage current of a corresponding clas- 
sical wave packet in the same billiard. In the classical 
calculations we take into account the diffractive effects 
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FIG. 3: Quantum mechanical transmitted current through a 
square billiard for different incoming modes of the wave packet 
a = 1, 2, 3 with (k)n/w = 1.5, 2.5, 3.5 respectively. The upper 
curve indicates corresponding classical current. Vertical dot- 
ted lines are guides for an eye indicating the same positions 
of the peaks in the current. Dashed lines give the asymp- 
totic power-law decay obtained by the best fit in the interval 
80 < t < 2000. The Heisenberg time th is indicated by a 
dot-dashed line. The curves are shifted for clarity. 



in the leads in the framework of the standard Fraunhofer 
diffraction approximation by injecting the electrons with 
the corresponding angular distribution P(9) (inset in Fig. 
U shows a calculated angular probability distribution for 
the lead geometry under consideration). A very good 
correspondence between the QM and the classical results 
allows us to ascribe each peak in the quantum-mechanical 
transmitted/reflected currents to a specific classical tra- 
jectory connecting the billiard leads, see inset in Fig. ||. 

A relative height of each peak depends on the density 
of the corresponding trajectories, and the angular distri- 
bution for a given incoming mode a. The effect of the 
later is clearly seen in Fig. |, where the QM current in 
a square billiard was shown for three different incoming 
modes of the wave packet a — 1,2, 3. At the initial stage 
of the current decay, t < 20, the positions of peaks are 
the same for all incoming modes, whereas their absolute 
values are different. This is explained by the fact that the 
angular distribution P(9) is different for different a, with 
its maximum being shifted to larger 9 for higher modes a. 
We conclude this discussion by noticing that all quantum 
billiards studied here exhibit similar characteristic peaks 
in the current at t < th that can be explained in terms 
of corresponding classical trajectories. 

Let us now focus on a long-time asymptotics of the 
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FIG. 4: Quantum mechanical reflected and transmitted cur- 
rent through a stadium-shaped billiard. Lower curves show 
corresponding classical currents. The wave packet enters the 
billiard in the second mode a = 2; (k) = 2.5w/ir. The inset 
show a long-time assymptotics of classical escape along with 
the example of a bouncing-ball trajectory. Dashed lines give 
the asymptotic power-law decay obtained by the best fit in 
the interval 50 < t < 2000. The Heisenberg time tb is indi- 
cated by a dot-dashed line. Time is measured in the units of 
the traversal time of the equivalent square of the same area. 
The curves are shifted for clarity. 



position, existence of periodic orbits, etc.). 

Let us now discuss wave-packet evolution in a stadium- 
shaped billiard. Classically, this billiard exhibits chaotic 
dynamics, and its classical escape rate shows fast expo- 
nential decay, see Fig. ^. The quantum mechanical de- 
cay of such a system is however, qualitatively different. 
It obeys a power law similar to the one observed for the 
square. The difference between the classical and QM de- 
cay asymptotics becomes clearly discernible at the time 
scale corresponding to ~ 50 — 100 classical bounces in a 
billiard. With this respect it is important to stress that 
the difference in the power-law and the exponential es- 
cape for classical regular and chaotic systems becomes 
discernible at a comparable time interval Note that 
the billiard at hand is designed in such a way that the 
classical escape through the left lead changes its asymp- 
totics from the exponential one to a slower power-law de- 
cay at t ~ 500. This behavior is caused by the bouncing- 
ball orbits which are accessible via the left lead only 
(see inset in Fig. |]). It is interesting to note that the cor- 
responding QM current through the left lead also starts 
to show slower decay at t ~ 500 in comparison to the 
right lead. We therefore speculate that, even though the 
long-time asymptotics of the QM and classical decays are 
qualitatively different, the QM decay still reflects some 
features of the underlying classical dynamics. 



wave packet dynamics. We start with a square billiard, 
see Fig. ||[ Its classical escape rate is independent of 
the number of modes in the leads and is well approx- 
imated by a power law ~ i~ 3 . On the contrary, the 
calculated quantum-mechanical decay, does depend on 
the number of modes in the leads and follows the power 
law decay with the exponents £ = 2.3,3,3.5 for modes 
a = 1,2,3 correspondingly. These values are somehow 
different from those expected from Eq. ([!]), £ = 2,3,4. 
(In a billiard with two leads a number of decay channels 
is given by M = 2a). 

One of the reasons for the above discrepancy may be 
related to the fact that in a billiard system the details 
of the coupling between the leads can be important for 
the selection of particular states that mediate transport 
through the system ]l7j (on the contrary, Eq. (|l|) corre- 
sponds to the case when all resonant states are excited 
with the same probability at t = ). Note that eq. (||) 
corresponds to the weak (tunneling) coupling between 
the leads and the dots in contrast to the regime of the 
open dot considered here. It is also worth to mention 
that Eq. (Q) is based on the random matrix theory and 
similar stochastic approaches ||, |l(], |n) . The predictions 
of these theories tend to be rather general in nature and 
they usually fail to account for specific features of the 
geometry under consideration (such as details of the lead 



The qualitative difference between the QM and clas- 
sical escape represents one of the main findings of the 
present work. We also find that the asymptotic decay of 
a quantum system obeys a power law that depends on 
the number of decay channels only, and is not sensitive 
to the nature of classical dynamics (chaotic or regular). 
This makes us conclude that quantum mechanics smears 
out the difference between classical chaotic and regular 
motion. With this respect it is important to stress that 
the difference in the classical decay rate in chaotic, regu- 
lar or mixed system is often used in various semiclassical 
approaches to describe observed transport properties of 
the quantum systems (statistics of the fluctuations, the 
shape of weak localization, fractal conductance fluctua- 
tions, etc.) [0 |6[ @, ||, |§] ■ We demonstrate however, that 
the crossover to the QM power law decay may occur at 
the same time scale when the difference between classi- 
cal regular and chaotic systems becomes discernible. Our 
findings thus strongly indicate that some of the SC pre- 
dictions should be used with certain caution or even be 
revised. Finally, the results reported in the present paper 
can be tested experimentally in the variety of systems in- 
cluding semiconductor quantum dots jjj , microwave cav- 
ities acoustical 20 and optical billiards 
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